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ABSTRACT 

We present the emission line profile models of hydrogen and helium based on the results 
from axisymmetric magnetohydrodynamics (MHD) simulations of the wind formed near the 
disk-magnetosphere boundary of classical T Tauri stars (CTTSs). We extend the previous out- 
flow models of 'the conical-shell wind' by Romanova et al. to include a well defined magne- 
tospheric accretion funnel flow which is essential for modelling the optical and near-infrared 
hydrogen and helium lines of CTTSs. The MHD model with an intermediate mass-accretion 
rate shows outflows in conical-shell shape with a half opening angle ~ 35°. The flow prop- 
erties such as the maximum outflow speed in the conical-shell wind, maximum inflow speed 
in the accretion funnel, mass-accretion and mass-loss rates are comparable to those found in 
a typical CTTS. The density, velocity and modified temperature from the MHD simulations 
are used in a separate radiative transfer model to predict the line profiles and test the con- 
sistency of the MHD models with observations. The line profiles are computed with various 
combinations of X-ray luminosities, temperatures of X-ray emitting plasma, and inclination 
angles. A rich diversity of line profile morphology is found, and many of the model profiles 
are very similar to those found in observations. We find that the conical-shell wind may con- 
tribute to the emission in some hydrogen lines (e.g. Ha, H/3, Pa/3 and Pa7) significantly when 
the temperature in the wind is relatively high (e.g. ~ 10 4 K); however, the wind contribution 
decreases rapidly when a lower wind temperature is adopted. The model well reproduces a 
relatively narrow and low-velocity blueshifted absorption component in He I A 10830, which 
are often seen in observations. 
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1 INTRODUCTION 

Mass outflows are commonly found during mass accreting phases 
of pre-main-sequence stars. While powerful collimated outflows 
are observed in protostars, often less collimated and weaker out- 
flows are associated with classical T Tauri stars (CTTSs) (e.g. Ray 
|et al.| [2007). The outflows in CTTSs are most likely powered by 
accretion processes, as the strengths of mass outflow indicators 
are known to correlate with those of mass accretion (e.g. |Cabrit| 
et al. 1990; Hartigan, Edwards, & Ghandour 1995). Understanding 
the origin and properties of the wind from young stellar objects is 
very important as it is closely related to their spin and disc evolu- 
tions. The mass-loss process is one of the key mechanisms through 
which YSOs can remove their angular momentum to the surround- 



ing environment during the active accretion phases (e.g. Hartmann 



& Stauffer 



2005,2007 



T989"l|K6nigl||199H |Shu et al.|[l994l |Matt & Pudritz 



2008). 
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The outflows in the low-mass stars such as CTTSs are un- 
likely driven by thermally or radiatively due to their relatively low 
temperature and luminosity. On other hand, relatively larg e mag- 



netic field strengths (~ 10 G) are found in CTTSs (e.g. Johns 



Krull et al.|1999l|Symington et al.|2003bl|Donati J.-R et al.|2007| 



Donati J.-F. et al.||2008| l which supports an idea that the wind is 
formed in magnetohydrodynamics processes (Konigl, Romanova, 
|& Lov elace 2011). The recent magnetohydrodynamics (MHD) 
simulation of Romanova et al. (2009) (hereafter RO09) showed a 
new type of outflow solution, so-called 'the conical-shell wind', 
which is produced when the stellar dipole magnetic field is com- 
pressed by the accretion disc into the X-wind (Shu et aTT||1994| > 
like configuration. The resulting outflow occurs in rather narrow 
conical-shell shapes with their half opening angles between 30° to 
40°. This conical-shell wind is similar to the X-wind of Sh u et al.l 
( 1994) in which the wind launching region is restricted to a narrow 
region near the inner edge of the accretion disc around the corota- 
tion radius. However, unlike the X-wind model, the conical-shell 
wind model of RO09 does not require the magnetospheric radius to 
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be matched with the corotation radius of the system. More recent 
simulations performed with larger radial domains have shown that 
the conical-shell winds become strongly collimated by the mag- 
netic force at larger distances (Lii, Romanova, & Lovelace|2012fr . 
The model has been also applied to the higher mass-accretion rate 
FU Orionis systems (FUORs) Konigl et al. (2011). Other possible 
wind configurations often considered for the outflows in CTTSs 
(e.g. Ferreira, Dougados, & Cabrit 2006) are: (1) a disc wind, 
which requires a sufficient magnetic field and ionization fraction 
to launch magnetocentrifugal winds over a relatively wide range 
of disk radii (e.g. Ustyugova et al. 1995; Romanova et al. 1997 



Ouyed & Pudritz 1997, Ustyugova et al. 1999; Konigl & Pudritz 



2000| IKrasnopol sky, Li, & Bland ford 2003] jPudritz et aL|20(77) , 
and (3) a stellar wind in which outflows occurs along the open 
magnetic fields from the stellar surface (e.g., Hartmann, Avrett, & 



Edwards 1982; Kwan & Tademaru 1 988]|Hirose et al.|1997||Ro 
manova et al. 2005 , Cranmer 2009). 

Spectroscopic studies of strong emission lines in CTTSs may 
be able to constrain the launching regions of the winds. This, in 
turn, would restrict possible wind formation theories mentioned 
above. A usefulness of the optically thick He I A10830 as a diag- 
nostic of the inner wind from the accreting stars has been demon- 
strated in the observations by e.g. Takami et al. (2002), Edwards 
|et al.| ( |2003),|Dupree et al.|f2 005 1, Edwards et al. (2006) (hereafter 
ED06) and Podio et al. (2008). In particular, ED06 showed that a 
very high fraction (~ 70 per cent) of CTTSs exhibit a blueshifted 
absorption (below the continuum), which is clear sign of outflow. 
This is much higher than the rate found in another wind sensitive 
line, Ha, for which only about 10 per cent of stars exhibit a sim- 
ilar type of blueshifted absorption component (below the contin- 
uum) (e.g., Reipurth, Pedrosa, & Lago 1996). Edwards (2007) and 
Kwan, Edwards, & Fischer (2007) suggested that the blueshifted 
absorption component in the He I A10830 profiles is cased by a 
stellar wind in about 40 per cent, and by a disc wind in about 
30 per cent of the samples in ED06. Using local excitation cal- 
culations combined with spectroscopic observations, |KwarT & Fis- 
|cher| ( |20iT) also demonstrated the effectiveness of He I A10830 for 
probing the density and temperature of the inner wind. The signs 
for the inner winds are also seen blueshifted absorption component 
in Na D, Ca II (H and K), and Mg II (h and k) (e.g. |Calvet| 1997) 
Ardila et al. 2002); however, probing the launching regions using 
these lines arc more difficult (ED06). 

Previously, we have developed and tested a multidimensional 
non-LTE radiative transfer model of hydrogen and helium line pro- 
files, which uses the Sobolev approximation in the source func- 
tion calculations (Kurosawa, Romanova, & Harries 2011, here- 
after KU11; see also Harries 2000, Kurosawa, Harries, & Syming- 
|ton|2006| >. The model uses rather simple kinematic and geometric 
models for both stellar and disc winds (see also Hartmann et al. 
1 1990| |Calvet, Hartmann, & Hewett||1992) . We demonstrated that 
our models are consistent with the scenario in which the narrow 
blueshifted absorption component of He I A10830 seen in observa- 
tions is caused by a disc wind, and the wider blueshifted absorp- 
tion component (the P-Cygni profile) is caused by a bipolar stel- 
lar wind. Similar conclusions were reached in the earlier study by 
|Kwan et aT| ( |2007| l who also applied simple kinematic and geomet- 
ric wind models. To advance our understanding of the line forma- 
tions in the inner wind of CTTSs, we now apply our model to more 
realistic wind configurations such as those found in MHD simu- 
lations. Therefore, our main goals in this paper are to present the 
optical and near-infrared hydrogen and helium line profile models 
computed with the density, velocity and modified temperature (Sec- 



tion BTTJ from the conical-shell wind solution found in the MHD 
simulations of RO09, and to check the consistency of the wind so- 
lution with observed line profiles. 

In Section [2] the model assumption and setup of our MHD 
simulations are presented, and those for the radiative transfer model 
for computing line profiles are given in Section [3] The results of 
the MHD simulations and the line profiles models are summarised 
in Sections H] and B] respectively. Brief discussions in which the 
line profile models are compared with observations are given in 
Section [6] Finally main findings and conclusions are summarised 
in Section|7] 



2 MHD MODEL 

The numerical method and the code used to simulate the out- 
flow from the disk-magnetosphere interaction region are essentially 
identical to those in the conical-shell wind models of Romanova] 
|et al.| p009l > (hereafter RO09; see also |Ustyugova eFaLl[2006J ); 
hence, we briefly describe only the most important aspects of our 
numerical MHD simulations below. 



2.1 Outline 

The simulation domain is divided into two regions: (1) the accre- 
tion disc region and (2) the corona region which are located outside 
of the disc region. The former is described by the equations of vis- 
cous and resistive MHD, and the latter by those of ideal MHD. 
The plasma is assumed to be an ideal gas with adiabatic index 



7 = 5/3. We adopt the standard a prescriptions of Shakura & Sun 



|yaev| (fl973) for the coefficient of the turbulent kinematic viscosity 
Ut = a v c s /O.K, where c s , Ok, a„ are the isothermal sound speed, 
the Keplerian angular velocity, and the dimensionless viscosity co- 
efficient, respectively. Similarly, the turbulent magnetic diffusivity 
coefficient is defined as ry t = adcf/fix where ay are the dimen- 
sionless turbulent diffusivity coefficient. Both, viscosity and diffu- 
sivity correspond to some type of turbulence in the disc; hence, they 
are included only in the disc which is distinguished from the low- 
density corona. Thus, we set a v , ad = when the density is lower 
than a typical density in the disc (pd), i.e. when p < pd- The vis- 
cosity term helps to bring matter towards the star in a steady rate, 
while diffusivity helps this matter to diffuse through the magnetic 
field lines of the magnetosphere. We refer to RO09 and |Lii et al.| 
(2012) for more comprehensive descriptions of our implementa- 
tions of the stress tensor, viscosity and diffusivity. 

The MHD equations are solved numerically by using the 
Godunov-type conservative code described in detail in RO09 and 
Ustyugova et al. (2006). The code uses the spherical coordinates 
(r, 6, (j>) where 9 and <j) are the polar and azimuthal angles. We 
consider axisymmetric approximation, i.e. d/d<f> = 0. 

2.2 Initial and boundary conditions 

The initial configuration of the magnetic field is set by a dipole field 
aligned with the stellar spin axis, i.e. 



B = [3(fJ.-r)r- fJ,r 2 ]/r 



(1) 



where fi is the magnetic dipole moment of the star. The simu- 
lation domain is initialised with a low-density, high-temperature 
plasma (corona) which is described as the following density and 
pressure distributions: p — p c exp[GAf*/(7?.T c r)] and p — 
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p c exp[GM t /(TZT c r)] where T c , p c and p c are the temperature, 
density and pressure at the outer boundary, respectively. The sym- 
bols G and M» are the gravitational constant and the stellar mass. 
The initial velocity of the corona is set to zero everywhere. The 
outer boundary r — Rout are divided into a disc region #d ^ ^ 
7r/2, and a corona region ^ 9 < 9 C \ with 9<± ~ 65°. At the begin- 
ning of the simulation (t — 0), there is no high density disc matter 
in the domain. We fix the density and temperature in the disc region 
at the outer boundary to p = pd and T = Td . Consequently as the 
simulation advances in time, a high-density (pd) low-temperature 
(Td) gas enters the simulation domain through the disc boundary 
region, 9 > 9^. The gas in the disc region continues to flow inward 
due to the a viscosity described earlier. The spin of the central star 
is initially set to a small value corresponding to r cor = -Rout where 
'"cor = (GM^/Q, 2 ) 1 ' 3 , and it is gradually increased to a final stel- 
lar angular velocity, fi, . As the simulation progresses, the informa- 
tion about the stellar rotation propagates rapidly to the low-density 
corona. 

At the inner boundary (r = Ri n ), the frozen-in condition 
is applied to the poloidal component B p of the field i.e., B r is 
fixed while 'free' boundary conditions are applied to Be and B,p 
(dBg/dr — and dB^/dr — 0). The free boundary conditions 
[d (■ • ■ ) /dr = 0] are also applied to density, pressure, and spe- 
cific entropy. The velocity components are calculated using free 
boundary conditions, and adjusted to be parallel to the magnetic 
field vector at the inner boundary. No matter enters from the inner 
boundary. The mass injection from the inner boundary, e.g. via a 
stellar wind, is not considered in this work. 

At the outer boundary (r = -Rout) of the corona region 
(0 ^ 9 < 9d), the free boundary conditions are applied to all 
the hydrodynamic variables. The matter is free to flow out, but no 
matter is allowed to enter from the outer boundary in the coronal 
region. At the outer boundary of the disc region, (0d ^ 9 ^ 7r/2), 
the density is fixed at the 'disc density' (pd). The velocity is fixed 
at a slightly sub-Keplerian velocity, i.e. ild = K,Q.(r c \) where 
ft = 0.997. This allows the matter to flow into the simulation region 
through the boundary more easily. Finally, the boundary conditions 
on the equatorial plane and on the rotation axis are symmetric and 
antisymmetric, respectively. 



3 RADIATIVE TRANSFER MODEL 

The density, velocity and modified temperature (Section [5.1^ from 
a snapshot of the MHD model, as described above, will be used 
as inputs for the line profile calculations. The hydrogen and he- 
lium line profiles from the simulations are computed by using the 
radiative transfer code TORUS (e.g. Harries 2000, Kurosawa et al. 
|2006[ KU11). The numerical method used in this work is essen- 
tially identical to those in KU11. Below, we briefly describe only 
the essential aspects of the model, and refer readers to KU1 1 for 
more comprehensive descriptions of our radiative transfer model. 



3.1 Outline 

The basic steps for computing the line profiles are as follows: 
(1) mapping of the density, velocity and temperature (but see Sec- 
tion ^. 1} , from an MHD simulation output to the radiative transfer 
grid, (2) the source function (S v ) calculation and (3) the observed 
flux/profile calculation. A brief description of each step is described 
below. 

In step (1), we use a mesh refinement grid which allows us an 



accurate mapping of an MHD simulation data on to the radiative 
transfer grid. In the grid construction process, we ensure that the 
resolution of the radiative transfer grid is higher or equivalent to 
the original MHD simulation grid. In step (2), we adopt the method 
described by Klein & Castor (1978) (see also Rybicki & Hummer 
|1978| |Hartmann, Hewett, & Calvet||1994) in which the Sobolev 
approximation method is used. Our atomic model consists of 20 
bound levels for H I, 19 levels for He I (up to the principal quantum 
number n = 4 level), 10 levels for Hell, and the continuum level 
for He III. In step (3), an observed line profile will be computed by 
using the source function found in step (2). The observed fluxes are 
computed in the cylindrical coordinates system whose symmetry 
axis points towards an observer. The line profiles are numerically 
computed via integrations of the equation of radiative transfer. The 
number of grid points used for a typical integration are n r — 180 
and n<t> — 100 for the radial and azimuthal grid points, respec- 
tively. A typical number of frequency points used across a given 
line profiles is n„ = 101. The line broadening effects (Stark and 
van der Waals) are also included for some hydrogen lines (e.g. Ha 
and H/3). 



3.2 Continuum sources 

3.2.1 Photosphere 

For all the line profile models presented in this work, we adopt 
stellar parameters of a typical classical T Tauri star for the cen- 
tral continuum source, i.e. its stellar radius R„ = 2.0 Rq and 
its mass M* = 0.8 Mq. Thus, we adopt the model atmosphere 
of Kurucz (1979), with the effective temperature of photosphere 
T p h = 4, 000 K and the surface gravity log g, = 3.5 (cgs), as the 
photospheric contribution to the continuum flux. 



3.2.2 Hotspots 

The hot spots on CTTSs are formed by the infalling gas along the 
magnetic field on to the stellar surface. As the gas approaches the 
surface, it decelerates in a strong shock, and is heated to ~ 10 6 K. 
The X-ray radiation produced in the shock will be absorbed by the 
gas locally, and re-emitted as optical and UV light (Calvet & Gull- 
bring 1998, Gullbring et al. 2000) - forming the high temperature 
regions on the stellar surface with which the magnetic field inter- 
sects. In this study, we adopt the hot spot temperature model de- 
scribed by Romanova et al. (2004) in which the local hot spot lu- 
minosity is computed directly from the inflowing mass flux on the 
surface of the star from MHD simulations. The temperature of the 
hot spots is determined by the conversion of kinetic and internal en- 
ergy of infalling plasma to a thermal radiation. Using the mass flux 
(p v r ) at a given location on the stellar surface, one finds the corre- 
sponding thermal radiation temperature as T ns = (pv r e 2 /a) 
where e g = (y 2 /2 + w) . Further, p, v r , v and a are the den- 
sity, the radial component of velocity, the speed of plasma, and the 
Stefan-Boltzmann constant, respectively. The specific enthalpy of 
the gas is k; = 7 (p/p) (7 — 1) where 7 is the adiabatic index and 
p is the gas pressure. We compare this temperature T ns with the 
effective temperature of photosphere T p h to determine the shape 
and the size of hot spots. When Th s > T pn , the location on the 
stellar surface is flagged as hot. For the hot surface, the total con- 
tinuum flux is the sum of the blackbody radiation with T ns and the 
flux from the model photosphere mentioned above. The contribu- 
tion from the inflow gas is ignored when Th s < T pn . 
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3.2.3 X-ray emission 

Relatively strong X-ray emission is commonly found in CTTSs, 
with their luminosities Lx (0.3 — lOkeV) ranging from ~ 
10 28 ergs _1 to ~ 10 31 ergs" 1 (e.g. |Telleschi et al.|2007][Gudel 



& Telleschi|2007||Gudel M. et al.|2010) . Recent studies by |Kwan 
& FischerH2011[ > and KU11 have shown that the photoionization 
by high energy photons, in addition to the normal continuum emis- 
sion from the photosphere, is most likely needed for generating 
the high opacities and hence strong absorption features seen in the 
observed He I A10830 line profiles of CTTSs (e.g. ED06). As in 
KU1 1, we adopt a simplest model of the X-ray radiation here. We 
assume that the X-ray radiation arises uniformly from the stellar 
surface as though it were formed in the chromosphere. This allows 
us to include the X-ray emission by simply adding it to the nor- 
mal stellar continuum flux (Section [3.2. 1) . The dependency of the 
line formation on the different X-ray emission locations is beyond 
the scope of this paper, but will be explored in the future. In this 
work, we assume that the X-ray emitting plasma is radiating ther- 
mally (as a blackbody) with a single temperature Tx. The thermal 
X-ray radiation flux is normalised with the total X-ray luminosity 
Lx (0.1 — 10 keV), which we set as an input parameter. Note that 
we restrict the contribution of the X-ray emitting gas to the contin- 
uum to be only within the X-ray energy range (0. 1-10 keV), i.e. the 
intensity (I„) from the X-ray emitting plasma is assumed to be zero 
for the energy below 0.1 keV and above 10 keV. The attenuation of 
X-ray flux is not taken into account in our model; however, this is 
to be improved in a future study. Further discussion on the X-ray 
optical depth can be found at the end of Section |5T2] 

Note that the stellar photospheric continuum in our model 
(Section |3.2. 1[ > has a relatively low effective temperature 
(4, 000 K); therefore, the stellar flux in EUV (including that at 
24.58 eV) is essentially zero. This is situation could be also thought 
as the case when the EUV flux is completely attenuated. With no 
EUV flux present, the only source of the photoionization of He I in 
our model is the X-ray emission as described above. 



4 MHD MODEL RESULTS 

4.1 Input parameters and grid setup 

The stellar parameters adopted here are M* = 0.8 Mq, R, — 
2.0 Rq and P — 5.38 d where M», R, and P are the stellar mass, 
radius and the rotational period, respectively. The corresponding 
corotation radius and the escape velocity of the system are 6.0 R t 
and 390 km s _1 , respectively. Following, RO09, we set the dimen- 
sionless viscosity and diffusivity coefficients (see Section [2~Tj > to 
q v = 0.3 and a<j = 0.1. For the models presented here, the surface 
magnetic field strength is fixed at _B* = 1 kG, and the magnitude 
of the dipole moment (equation [TJ at fj, — 40 in the units of B, . 
The initial corona temperature is set as T c = 9.2 x 10 6 K, and the 
disc temperature at the outer boundary is fixed at T<j = 1.5 x 10 3 K 
(see Section [2T) . The ratio of the corona to disc densities (p c /p c ) 
are fixed at ~ 1.7 x 10~ 4 , and we vary pd to control the disc mass- 
accretion rate. 

The simulations are performed in the radial range, R[ n ^ r ^ 
Rout where R[ n = 1.0 R* and Rout — 33.7 Ri n , and in the po- 
lar angle range, ^ 9 ^ 7r/2. The grid spacing in 9 is uniform, 
but that in r is chosen such that each side of the curvilinear rect- 
angle cells is approximately equal. The number of grid points used 
in the models are ng = 40 and n r — 96. Note that the grid setup 
used here is slightly different from those in the conical-shell wind 



models of RO09 who used the radial range: Ri n — 2.0 R, and 
Rout — 16 7?in, and the number of grid points: ng = 31 and 
n r — 51. Their corotation radius for a CTTS model is 3.0 R* 
which is 1/2 of our value. We have run simulations with higher 
resolutions (e.g. ng = 80 and n r = 192), but overall flow patterns 
as well as the physical quantities (e.g. plasma temperature, density 
and velocity) are very similar to those of the lower resolution mod- 
els presented in this work. 

4.2 Matter flow in conical-shell winds 

We have performed the MHD simulations with ranges of the disc 
density pd and have examined the dependency of the flows on 
Pd. Here, we present two most representative cases for (1) a rel- 
atively high mass-accretion rate (Model A) and (2) a moderate 
mass-accretion rate (Model B). Fig. [TJ shows the maps of den- 
sity (p), poloidal mass-flux ($ m = pv p ) along with sample mag- 
netic field lines and the directions of the poloidal velocity (ti p ), for 
both models. The main model parameters and the corresponding 
mass-accretion rates (M a ), mass-loss rates in the conical-shell wind 
(Mew) and the mass-loss rates in the polar wind (M pw ) are sum- 
marised in Tab. [TJ In the following, each model will be described 
in more detail. 

4.2.1 A relatively high mass-accretion rate case: Model A 



This model uses a relatively large disc density (at the outer bound- 
ary), pd = 4.9 x 10~ 12 gcm -3 . The disc matter initially en- 
tering from the outer boundary forms an accretion disc around 
the star, and reaches a semi-steady state as shown in Fig. [Tj The 
high disc density matter pushes the stellar dipole magnetic fields 
almost all the way to the stellar surface, and a very small fun- 
nel (~ 1.2 7?* in size) is formed. The matter accretes on to the 
stellar surface through the small funnel with the accretion rate 
M a = 3.4 x 10" 7 MQyr _1 . The geometrically inclined field 
lines near r ~ 2 R, on the disc plane and the inflated field above 
it create conditions favourable for matter from the inner disc re- 
gion (RO09). The gas flows, from the inner disc region, into a 
conical-shell shaped wind. The flow is somewhat collimated near 
the outer boundary and has a half-opening angle 9 ~ 18°. This 
is smaller than those of the conical-shell wind simulations of RO09 
(6* ~ 30° — 40°) who used a smaller radial range (Rout = 16 Ri n ) 
than our value (Rout = 33.7-Rin). On the other hand, [Lil et al.| 
(2012) obtained a much more collimated outflow (9 ~ 4°) using 
a larger radial range (Rout = 42 Ri n ). 

The poloidal and toroidal velocity components (v p and v^) 
along with the magnitude of the total velocity (t>tot) along the 
conical-shell wind (along the black dashed line in the upper right 
panel of Fig.[TJ( are shown in Fig. [2] The rotational speed v$ domi- 
nates near the base of the wind, and its value is essentially the same 
as that of the Keplerian velocity at the wind base. The poloidal 
component v p increases along the conical-shell wind, and it peaks 
around d ~ 8 R, with v p ~ 125 km s" 1 , then it slowly decreases 
beyond that point, and v p reaches ~ 60 km s _1 when the conical- 
shell wind escapes from the outer boundary. The slight decrease in 
tip at large radii is related to the fact that the reference lines in Fig.[Tj 
are not exactly parallel to the outflow stream lines, and v p tends to 
decrease as the distance from the symmetry axis (z) increases at a 
given value of z. The velocity range found here is consistent with 
the extent of the narrow blueshifted absorption components (up to 
~ 200 km s _1 ) in the observed He I A10830 profiles obtained by 
e.g. ED06. 
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Table 1. Summary of MHD Models 



Line profile models of the inner wind ofCTTSs 



Model ID 



'id 
(K) 



(K) 



Pd 

(gem -3 ) 



Pc 



ilia 



Al c 



My, 



(gem" 3 ) (Mgyr 1 ) (Mgyr 1 ) (Mgyr" 1 ) 



Funnel Flow Size 



A 1.5 x 10 3 9.2 x 10 6 4.9 x 10" 12 8.2 x 10~ 16 
B 1.5 x 10 3 9.2 x 10 6 8.2 x 1CT 13 1.4 x 1CT 16 



3.4 x 1CT 7 
4.1 x 10~ s 



8.2 x 1CT 8 
6.9 x 10~ 9 



1.4 x icr 10 
2.6 x icr 10 



< 1.2 
-2.0 






Figure 1. Summary of MHD models with a high mass-accretion rate (Models A; top panels) and a moderate mass-accretion rate (Model B; bottom panels). 
See Table[T]for the model parameters. The density maps (in logarithmic scale; in cgs units) are over-plotted with the directions of the poloidal velocity (v p ) as 
arrows (left panels). The poloidal matter flux v p p maps (in logarithmic scale; in cgs units) are shown with sample magnetic field lines (middle panels). The 
zoom-in poloidal matter flux v p p maps (right panels) are shown with both sample magnetic field lines and poloidal velocity directions as arrows. In the left 
panels, the reference lines passing through the conical-shell (black dashed line) and polar winds (green dashed line) are also shown. The length scales are in 
the units of the stellar radius (H*). 



Fig. [T] also shows the lower density and higher speed out- 
flow in the polar direction which we refer to as 'the polar wind.' 
This should be distinguished from 'the stellar wind' which arises 
from/near the stellar surface, but is not implemented in this model. 
The matter in the polar wind is simply the gas redirected from the 
upper part of the accretion funnel flow, and it does not originate 
from the stellar surface. The poloidal component of the velocity 
in the polar wind (along the green dashed line in the upper right 
panel of Fig. fTl is also shown in Fig. [2] The polar wind quickly 
accelerates to about 1.3 x 10 3 km s~ in a few stellar radii, and 
it slowly decelerates beyond d ~ 15 R*. The density in the po- 
lar wind is very low and the polar region is magnetically domi- 
nated. As a result, the low density gas is accelerated up to a very 
high velocity via strong magnetic forces. A similar component of 
the outflow has been observed in the simulations for the propeller 
regime (a fast stellar rotation regime) (e.g. Ustyugova et al. 2006, 
|Romanova et al.|2 009). As mentioned before, the slight decrease in 



v p at large radii seen in Fig. [2] is related to the fact that the refer- 
ence lines in Fig. [T] are not exactly parallel to the outflow stream 
lines. This velocity is much larger than the extent of the wide and 
deep blueshifted absorption feature (~ 400 km s _1 ) in the ob- 
served He I A10830 profiles (e.g. ED06). The poloidal velocity of 
the polar wind is about 10 times larger than that of the conical- 
shell wind in most of the radial ranges. The mass-loss rates in the 
conical-shell and polar winds are A/ cw = 8.2 x 10 _8 Mq yr -1 
and M pw = 1.4 x 10 _10 Mq vr" 1 , respectively. The former is 
much larger and the latter is smaller than a typical mass-loss rate of 
CTTSs (~ 10~ 9 Mq yr -1 ), but they are still within the range of 
the observed mass-loss rates (e.g. Hartigan et al. 1995). 

Although the mass-accretion rate M a = 3.4 x 10 -7 Mq yr -1 
found in Model A is much higher than that of a typical CTTS 



10~ 8 Mq yr _1 (e.g. Gullbring et al 



1998,Hartmannetal. 



1998 



|Calvet, Hartmann, & Strom 2000), it is still an acceptable value 



for CTTSs (e.g. Hartigan et al. 1995). However, this model is not 
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Figure 2. The total speed (ftot; dash-dot), the poloidal (v p ; solid) and azimuthal components (y$; dash) of the velocity along the reference lines indicated in 
the left panels of Fig.fTlare shown for Models A (top panels) and B (bottom panels). The horizontal axis indicates the distance (d) from the base of the wind 
along the reference lines. The panels on the left column are the velocity components along the reference line that goes through the 'conical-shell winds.' The 
panels on the right column show v p along the reference line that go through 'the polar wind'. Compared to v p , v^, is negligibly small and v p rs utot in the 
polar wind; hence, they are omitted from the plots. 



suitable for modelling line profile for CTTSs because the size of 
the magnetospheric accretion funnel flow is much smaller than 
that indicated by previous studies (e.g. Muzerolle, Calvet, & Hart- 
|mann|1998a||Muzerolle, Calvet, & Hartmann 200J}. The observed 
line profiles (e.g. Pa/3, Br7 and He I A10830) often show red- 
shifted absorption features that often extends to 200-300 km s _1 
(e.g. |Edwards et al.|1994l|Alencar & Basri|2000| |Folha & Emer-| 
son 2001 1 which are formed in the inflowing gas that is accreting 
on to the stellar surface through the funnel flows. When the size 
of the magnetosphere or the funnel flow is too small, the inflow 
speed in the funnel also becomes too small since the gas gains the 
kinetic energy from the gravitational potential energy (c.f. Ghosh, 
|Pethick, & Lambl[T977l |Hartmann et~aLl[T994] . In Model A, we 
find that the poloidal velocity of the gas in the funnel flow is only 
~ 120 km s _1 which is much smaller than a typical observed value 
(200-300 km s _1 ). Further, the geometry of the funnel flow is not 
ideal for producing the redshifted absorption feature. Since the fun- 
nel flow is located very close to the stellar surface, the funnel stream 
is almost parallel to the stellar surface. This flow geometry reduces 
the line-of-sight (recession) speed of the gas toward an observer. 
The line profile models based on this type of flow will result in a 
very weak redshifted absorption component with a velocity extent 



that is too small when compared to observations. To improve the 
MHD model to be more suitable for CTTSs, we now attempt to 
enlarge the size of the accretion funnel by preventing the accretion 
disc from pushing the stellar dipole magnetic field too close to star. 



4.2.2 A moderate mass-accretion rate case: Model B 

To reduce the amount of the compression of the magnetosphere to- 
wards the stellar surface by the accretion disc, the disc density pa 
at outer boundary is now reduced by a factor of ~ 6 from the value 
used in Model A. The lower density plasma reduces the total gravi- 
tational force on the magnetic field lines, and consequently the size 
of the magnetosphere or the accretion funnel becomes much larger 
(~ 2.0 R,) than the pervious case, Model A. In Model B, the fun- 
nel flow is clearly detached from the stellar surface (Fig. fTl. The 
corresponding mass accretion rate on to the stellar surface through 
the funnel is M a = 4.1 x 10~ 8 Mq yr _1 , which is very sim ilar to 
the value found for a typical CTTS (~ 10" 8 M yr" 1 , e.g.|Gull 



|bring e t al. 1998] |Hartmann et al.p998J|Calvet et~aT]|2000f , and 
is about ~ 8 times smaller than that of Model A. The inner edge 
of the accretion disc is now located ~ 2.0 R, from the origin, and 
the outflow in the conical-shell wind originates at a larger distance 
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r w 4.0 R, (in the disc plane) which is slightly smaller than the 
corotation radius (R cr — 6.0 i?*). 

The outflow geometry is similar to that of Model A, but the 
wind is less collimated with its half-opening angle of the conical- 
shell wind (measured at the outer boundary) O ~ 35° (Fig. [TJ, 
which is similar to the value found in RO09. Fig. [2] also shows 
the poloidal and toroidal velocity components (v p and v$), and the 
magnitude of the total velocity (vtot) along the conical-shell wind 
(along the black dashed line in the lower right panel of Fig.fTl. As 
also seen in Model A, the rotational speed V4, dominates near the 
base of the wind. The poloidal component v p increases along the 
conical-shell wind, and it peaks around d (a 15 Rq with v v ~ 
125 km s" 1 . The flow then slowly decelerates beyond that point, 
and v p reaches ~ 90 km s" 1 near the outer boundary. As in the 
previous model (Model A), the slight decrease in v p at large radii 
is related to the fact that the reference lines in Fig. [Hare not exactly 
parallel to the outflow stream lines. 

Again, this velocity range in the conical-shell wind is con- 
sistent with the velocity range of the narrow blueshifted absorp- 
tion components (up to ~ 200 km s" 1 ) found in the observed 
He I A10830 profiles (e.g. ED06). Note that the outflow speed (v p ) 
in the conical-shell wind in Model B is similar to that of Model 
A. The polar wind quickly accelerates to ~ 550 km s" 1 (along the 
green dashed line in the lower right panel of Fig.fTl at the distance 
d « 5 R* from the base of the wind. The wind speed then slightly 
decreases, but it remains almost constant all the way to the outer 
boundary. The outflow speed (v p ) of the polar wind in Model B is 
about 2 times less than that of Model A. This velocity is compara- 
ble to the extent of the wide and deep blueshifted absorption feature 
(~ 400 km s" 1 ) in the observed He I A10830 profiles (e.g. ED06). 
The mass-loss rates in the conical-shell and polar winds are M cw = 
6.9 x 1O" 9 M yr" 1 and M pw = 2.6 x 10~ 10 Mq yr~\ respec- 
tively. The former is larger and the latter is smaller than a typi- 
cal mass-loss rate of CTTSs (~ 10" 9 M yr" 1 ), but they are still 
within the range of the observed mass-loss rates (e.g. Hartigan et al. 
[T995) . 

The geometry of the well defined accretion funnel in Model B 
is very similar to that of the magnetospheric accretion funnels used 
in hydrogen line profiles models e.g. Hartmann et al. ( 19941, Muze- 
rolle et al. (20011, Symington, Harries, & Kurosawa (2005a) and 



Table 2. Summary of line profile model parameters 



Kurosawa et al. 2006) who adopted the axisymmetric magneto 
sphere model of Ghosh et al. ( 19771 (see also |Ghosh & Lamb| 
1979a. Ghosh & Lamb 1979b I. These models have been success- 
fully demonstrated that the redshifted absorption features or the in- 
verse P-Cygni profiles seen in some of the observed hydrogen lines 
can be explained by the absorption by the inflowing gas through the 
magnetospheric accretion funnels. In Model B, we find that the flow 
speed in the funnel near the stellar surface reaches ~ 250 km s - , 
which is very sim ilar to the values seen in the observation, 200- 
300 kms -1 (e.g. Edwards et al.||l994 Muzerolle, Hartmann, & 



|Calvet|1998bl|Alencar & Basri|2 000. Folha & Emerson|2001 

In summary, Model B shows very similar inflow and outflow 
velocities, mass-loss and mass-accretion rates to those found in the 
observations. Hence, it is a very good candidate for being tested 
with our radiative transfer model which uses the MHD simulation 
results as model inputs. In the following, we describe the hydrogen 
and helium line profile models based on the flow solution found in 
Model B. 
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5 LINE PROFILE MODELS 

We now investigate whether the MHD outflows found in Model B 
of the previous section is a plausible solution for CTTSs, by com- 
puting line profiles based on the simulation, and by checking the 
consistency with the existing spectroscopic observations. The stel- 
lar parameters adopted here are same as in Section |4.1| and the 
continuum sources as described in Section [3~2| The density, ve- 
locity and temperature (but see Section [5~T) found in Model B are 
mapped on to the radiative transfer grid (e.g. KU1 1), and used as in- 
put variables in the profile calculation. We exclude the high density 
disc accretion disc region (indicated by the dashed line in Fig. 13), 
and concentrate on modelling the line emission and absorption only 
from the accretion funnel and the outflow regions. The disc region 
is excluded from the computations because of the gas temperature 
near the inner part of the disc, especially near the inner edge of the 
disc, is too high (T ~ 10° K) in the simulation result from Model B 
(Fig. 13). If the relatively high temperature gas near the inner edge of 
the disc is included, the line emission is dominated from this region, 
and the line strengths are unrealistically strong. The disc tempera- 
ture at larger radii (r > 20 R t ) are low 1,500-2,000 K, and we do 
not expect any significant contribution to the atomic hydrogen and 
helium emission. To mimic the opaqueness of the accretion disc, 
we simply place a geometrically thin but optically thick disc on the 
equatorial plane with the inner radius corresponding to the inner ra- 
dius of the high density MHD disc region (the vertical dashed line 
in Fig. [3) which is removed from the radiative transfer calculation. 
In the following, we examine the dependency of the line pro- 
files on our main model parameters: the maximum gas temperature 
(r max ), X-ray luminosity (Lx) and the temperature of the X-ray 
emitting gas (Tx) (cf. Section [3. 2. 3\ . The model parameters are 
summarised in Table [2] To examine the dependency of the line 
profiles on the model parameters, we mainly focus on modelling 
the lines which often show a sign of wind (a blueshifted absorp- 
tion component), namely Hq and He I A10830. In the following 
two sections (Sections |5 . 1 1 and |5.2[ l, example line profiles are cal- 
culated at i = 30° because He I A10830 profile exhibits both the 
blueshifted wind absorption and the redshifted magnetospheric ab- 
sorption most clearly at this inclination angle. However, readers 
should be aware that the probabilities of observing an object with 
i < 30° and i > 30° are 0.134 and 0.866, respectively. Exam- 
ple calculations of additional hydrogen and helium lines, and their 
dependency on the inclination angle (i) will be also presented in 
Sectionl5.3l 



5.1 Dependency on the maximum gas temperature 

Our initial calculations of the line profiles based on the MHD sim- 
ulation Model B, with unmodified temperature structure from the 
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Figure 3. The temperature maps (in logarithmic scale; in K) of MHD Model B (cf. TablefTland Fig.fTJ in full-scale (left panels) and in small-scale {right 
panels). The upper panels show the original temperatures (unmodified) from the MHD simulation. The lower panels show the temperatures from the same 
MHD simulation, but the upper limit of the temperature is set to T max = 10 4 K. N.B. the different temperature ranges are used in the colour indeces for the 
upper and lower panels. The temperature maps from the lower panels are adopted in the line profile calculation (see Section [5TT) . The high density disc zone, 
indicate by the dashed lines is excluded from our line profile calculations. See SectionlTTlfor explanations. The length scales are in the units of the stellar 
radius (R*). 



simulations, showed that the line strengths (e.g. Ha) are unrealis- 
tically too strong when compared with observations (e.g. Reipurth 
|et al. 1 1 9961. This is mainly due to the emission from the relatively 
high temperature gas (~ 10 5 K) near the interface between the rel- 
atively high density conical-shell wind region and the low density 
polar wind region (Fig.[3](. The original temperature of the flows in 
the outflows are, in general, higher than the typical wind temper- 
atures used in our previous line profile models (~ 10 4 K, KU11) 
which showed reasonable agreement with the observed CTTSs. A 
recent local excitation study for line opacities of CTTSs by [Kwan| 
& Fischer (201 1 ) has also found a similar temperature range. When 
the gas from the inner part of the disc is ejected to the wind, it cools 
down due to the adiabatic expansion, which is treated properly in 
the MHD simulations. However, the additional radiative cooling 
mechanism, which is not included in the simulation, may be re- 
quired to cool down the gas even further. In general, the temper- 
atures from the MHD simulations could be quite uncertain. We 
therefore modify the original gas temperature of the MHD simu- 
lations, in order to match observed line profiles, and treat the tem- 
perature as though a free parameter in the line profile calculations. 
For this reason, we limit the maximum temperature (T max ) of the 



gas from the MHD simulation to be ~ 10 4 K when it is applied to 
the radiative transfer models. 

The lower panels in Fig. [3] show the temperature maps of 
MHD Model B when the maximum allowed temperature T max = 
10 4 K is applied. The figures show that T max not only affects the 
low density and high temperature polar wind regions, but also affect 
the conical-shell wind region. The relatively high density conical- 
shell wind region (including its base), is essentially isothermal 
(with T = T max = 10 4 K) in this case. On the other hand, the 
original temperature of the accretion funnel from the MHD simu- 
lation (the upper panels in Fig.[3](, is relatively low compared to the 
wind temperature and is mostly less than ~ 10 4 K. However, the 
temperatures in the funnel flow near the stellar surface and near the 
accretion disk are slightly higher than 10 4 K (in the original tem- 
perature from the MHD simulation); hence, they are affected by 
the imposition of Tmax. In the following, we examine the effect of 
changing T max on the computed line profiles. 

Fig.|6]shows the model line profiles for Ha and He I A10830 
computed for T max = 8.0 x 10 3 , 1.0 x 10 4 and 1.4 x 10 4 K 
(Models Bl, B2 and B3 in Table 12). The inclination angle, the 
X-ray luminosity and temperature are fixed at i — 30°, Lx — 
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4 x 10 30 erg s 1 and Tx = 2 x 10 6 K, respectively. The strengths 
of the emission lines are comparable to those found in observations 



(e.g. |Reipurth et al.||1996[|Edwards et al.||1994||Muzerolle et al.| 



|1998b||Alencar & Basri|2000| ED06) in this temperature range. 

The line profile morphology of the He I A10830 models 
are also similar to some of the observed line profiles found in 
ED06. In particular, they are similar to those which exhibit a 
narrow blueshifted absorption component that often extends to 
~ 200 km s . The relatively narrow blueshifted absorption com- 
ponents, centred around v w —120 km s , seen the models are 
caused by the outflowing gas in the conical-shell wind, and the red- 
shifted absorption component, which shows the minimum flux at 
V ~ 100 km s -1 and extends to v ~ 250 km s -1 is caused by 
the inflowing gas that is accreting on to the stellar surface through 
the accreting funnel. The figure shows that the emission strength 
and the absorption components (in both red and blue wing) of 
He I A 10830 are rather insensitive to the gas temperatures in this 
range. This is because of the relatively large ionization potential 
(24.6 eV) of He I, and the plasma temperature is not quite high 
enough for significant collisional excitations. The main cause of 
the excitation for He I in these models are due to photoionizations 
e.g. by the X-ray flux (c.f. KU1 1). 

On the other hand, the emission strengths of the model Ha 
is sensitive to T max value. The line emission becomes stronger 
as T max increases. Relatively weak but clearly visible absorp- 
tion components are seen in the blue wing of the line (v ~ 
— 120 km s -1 ). This, again, is caused by outflowing in the conical- 
shell wind. Unlike He I A10830, the redshifted absorption compo- 
nent that goes below the continuum level is absent from the Ha 
profiles. This is consistent with observations (e.g. Reipur ffTet al.| 
1996) which find that Ha profiles of CTTSs rarely show redshifted 
absorption below the continuum. For T max — lAx 10 4 K, the line 
centre flux reaches ~ 30 (normalised to the continuum) which is 
close to the maximum value found in the observations. The peak 
flux of the Ha profiles computed with T max = 1.0 x 10 4 K 
is ~ 15 (normalised to the continuum) which is similar to typi- 
cal values found in observations. In summary, based on the line 
strengths of Ha in Fig. El in comparison with those of typically 
observed Ha, the gas temperature of the accretion funnel and the 
conical-shell wind in this particular MHD simulation should be 
less than ~ 1.4 x 10 4 K. In the rest of this paper, we adopt 
Tmax = 1.0 x 10 4 K in our line profile calculations. 

Here, we have adopted a rather simplistic method of the gas 
temperature assignment, with the introduction of T max parameter. 
However, this leads to a large reduction of temperatures (from ~ 
10 6 K to ~ 10 4 K) in the low-density polar wind region where 
the gas temperature may be less affected by the radiative cooling. 
The limiting the gas temperature to ~ 10 4 K in the polar wind 
region would be reasonable when the wind density there is very 
low; hence, the omission of the opacity in the polar wind region 
does not affect the line profiles. On the other hand, if the polar wind 
density becomes relatively high, the wind opacity may become high 
enough to affect the line profile; therefore, the temperature of the 
polar wind should be treated more properly. 

5.2 Dependency on X-ray luminosity and temperature 

Next, we examine the dependency of the model profiles on the X- 
ray luminosity (Lx) and the temperature of the X-ray emitting gas 
(Tx). As in Section [3.2.3| our models assume that the X-ray emit- 
ting plasma is radiating thermally (as a blackbody) with a single 
temperature Tx. The thermal X-ray radiation flux is normalised 
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Figure 4. Dependency of He I A 10830 (left) and Ha (right) on the max- 
imum gas temperature allowed (T max ) in the profile calculations. MHD 
Model B is used for all the models shown here. The profiles are computed 
at the fixed inclination angle i = 30° and fixed X-ray luminosity and tem- 
perature, which are Lx = 4 x 10 30 ergs" 1 and Tx = 2 x 10 6 K, respec- 
tively. The models shown here use T max = 8.0 X 10 3 (solid), 1.0 X 10 4 
(dash), and 1.4 X 10 4 K (dash-dot), which corresponds to Models Bl, B2 
and B3 in Tablep] respectively. 



with the total X-ray luminosity Lx (0.1 — 10 keV). The photoion- 
ization process by high energy photons are expected to be important 
for the formation of the wind sensitive line, He I A10830 (c.f. |Kwaii| 
|&Fischer|20"TT] KUll). 

Fig. B] shows the model He I A10830 profiles computed for 
L x =4 x 10 29 , 4 x 10 30 and 2 x 10 31 ergs -1 (Models B4, B2 and 
B5 in Table [21, which are within the range of observed values for 
CTTSs (e.g. Tell eschi et al.|2007[|Gud"el & TelleschipOOTl |Gudel| 
|M. et al.||2010| l. The inclination angle and the X-ray temperature 
are fixed at i = 30° and Tx — 2 x 10 6 K, respectively. The corre- 
sponding photoionization rates of He I from the ground state (7n e i) 
at a typical location in the wind (r = 5.6-R,) for Models B4, B2 
and B5 are 1.2 x 10" 5 , 1.2 x 10" 4 and 6.0 x 10" 4 s -1 , respec- 
tively (Table [51. Since the X-ray flux is directly proportional to a 
value of Lx, the photoionization rate 7h c i is also proportional to 
Lx (see Table[2|. The photoioniza tion rates here are similar to those 



found in 



Kwan & Fischer 



2011i, i.e. 7hc 



10" 



10" 



the local excitation calculation of the plasma, in the context of the 
classical T Tauri stellar wind. 

The X-ray luminosity between 4 x 10 29 and 2 x 10 31 ergs -1 
(or equivalently 7 Hc i between 1.2 x 10 -5 and 6.0 x 10 -4 s -1 ) 
produces resonable line strengths in He I A10830 when compared 
to observations (e.g. ED06), and the line is very sensitive to this 
parameter. As in the previous tests for T max (Fig. HI, the model 
He I A10830 profiles show both the redshifted absorbtion casued 
by the accretion funnel and the blueshifted absorption component 
casued by the conical-shell wind. The figure shows that the line 
centre flux becomes larger as Lx (or equivalently Thsi) increases. 
Similarly, the depths of the absorption components on the blue 
sides become deeper and wider as the value of Lx (or equiva- 
lently 7hci) increases. Interestingly, for the model with the largest 
Lx (2 x 10 31 ergs -1 ), the blue absorption component extends to 
~ —200 km s -1 which is larger than a typical maximum speed ob- 
served in the conical-shell wind (see Fig.[2l. Although the deepest 
and main part of the blueshifted absorption (~ —120 km s -1 ) is 
casued by the relatively high density conical-shell wind, the shal- 
lower but higher velocity absorption is cased by the lower density 
and higher velocity (Fig. [2l gas in the polar wind. For the largest 
Lx model, there are enough high energy photons to photoionize 
He I in the polar wind, and the He I A10830 opaicty in the polar 
wind becomes non-negligible. 

The figure also shows the model Ha profiles computed for the 
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same set of Lx values. Unlike He I A 10830, Ha is insentitive to 
the value of Lx, and the computed profiles for all three Lx val- 
ues are almost indentical to each other. The model Ha profile does 
not depend on the value of Lx mainly because the photoioniza- 
tion rates for 'H I' in the X-ray energy range (0.1 — 10 keV) used 
here is small compared to the collinsional rates in the relatively 
high density gas in the funnel flow. A weak blueshifted absorption 
component, casued by the conical-shell wind, is also seen at v ~ 
— 120 km s _1 . Unlike the blueshifted absorption component seen 
in the He I A 10830 model with the largest Lx, the wind absorp- 
tion here does not extend to a high velocity (e.g. ~ —200 km s _1 ); 
hence, the polar wind is not contributing to the absorption in the 
Ha model profiles here. 

We now examine the effect of changing the X-ray temperature 
Tx on the profiles. Fig. [6] shows the model He I A 10830 profiles 
computed for Tx =2 x 10 6 , 5 x 10 B and 2 x 10 7 K (Models B2, 
B6 and B7, respectively in TablefTl. The inclination angle and the 
X-ray luminosity are fixed at i — 30° and Lx — 4 x 10 30 ergs -1 
for all three models. The corresponding photoionization rates of 
He I (7h c i) at a typical location in the wind (r — 5.6 7?*) for Mod- 
els B2, B6 and B7 are 1.2 x 10~ 4 , 2.0 x 10~ 5 and 1.2 x 10~ 6 s _1 , 
respectively (see Table [2l. For a fixed Lx value and in this tem- 
perature range, the photoionization rate of He I (7H01) decreases 
as Tx increases because the peak of the blackbody radiation curve 
shifts to a higher wavelength. Since the X-ray flux is normalized 
to a fixed value of Lx, this leads to a reduction of the softer X-ray 
flux when Tx is increased. Consequently, the photoionization rate 
of He I decreases because its photoionization cross section is much 
larger towards the EUV frequency. 

The line strengths of the model He I A10830 profiles computed 
with this range of Tx are comparable similar to those found in the 
observations (e.g. ED06). The line is sensitive to the change in the 
value of Tx (or equivalently 7Hei). The flux at the line centre de- 
creases as Tx increases (or equivalently as 7h c i decreases), and 
both blueshifted and redshifted absorption components becomes 
weaker as Tx increases (or equivalently 7H01 decreases). Similar 
to the earlier cases (Fig.B), the model Ha profiles are not sensitive 
to the change in Tx value, i.e. the model profiles are almost iden- 
tical to each other. Again, the line is insensitive to the change in 
Tx or more precisely to change in the X-ray flux because the pho- 
toionization rates for H I in the X-ray energy range (0.1 — 10 keV) 
is small compared to the collisional rates. 

As mentioned earlier (Section |3.2.3| l, the attenuation of the X- 
ray flux is neglected in our level population calculations; hence, it 
is implicitly assumed that the X-ray radiation is optically thin. To 
check the validity of this assumption, we have calculated the X- 
ray optical depth (rx) at the photon energies hv = 0.1, 1.0 and 
10.0 keV along the line of sight to an observer located at the incli- 
nation i = 60°. We found r x « 90, 0.4 and 0.01 for hv = 0.1, 1.0 
and 10.0 keV, respectively. This suggests that the X-ray flux should 
be almost completely attenuated for all photons with hv < 1.0 keV. 
On the other hand, the X-ray flux at hv > 1.0 keV are unaffected 
by the X-ray opacity. 

To assess the effect of the X-ray attenuation on the model pro- 
files, we do the following: (1) compute the continuum optical depth 
rx at all frequency in X-ray, (2) attenuate the X-ray flux by a fac- 
tor of e _TX at each X-ray frequency, (3) recompute line profiles 
(including source function calculations) using the X-ray flux cor- 
rected for the attenuation due to X-ray opacity, and (4) compare the 
new line profiles with the models computed without the attenua- 
tion of the X-ray flux. For this test, we use He I A10830 line only 
since Ha profiles are rather unaffected by the X-ray flux as previ- 



ously seen in Figs. Bland [6] We have recomputed the He I A10830 
line profiles for Models B2 and B6 (Table |2l with the X-ray flux 
with the attenuation. Note that these two models have the X-ray 
luminosity Lx = 4 x 10 30 ergs _1 but have different X-ray tem- 
peratures, Tx = 2 x 10 6 and and 5 x 10 6 K, respectively. The 
results are shown in Fig. [7] 

The figure shows that the reduction of the line emission and 
absorption strengths due to the attenuation of the X-ray flux is 
significant, for a relatively low X-ray temperature (Model B2: 
Tx = 2 X 10 6 K). However, the effect of the X-ray attenuation 
becomes notably smaller for the model with a higher X-ray tem- 
perature (Model B6: Tx = 5 x 10 6 K). The difference is mainly 
caused by the difference in the flux distributions of the blackbody 
radiation curve for different X-ray temperatures. The X-ray flux 
of the lower X-ray temperature (Tx = 2 x 10 6 K) peaks at a 
lower energy {hv ~ 0.5 keV) than that of the higher temperature 
(Tx =5x 10 6 K) which peaks at hv » 1.2 keV. As mentioned 
earlier, the X-ray attenuation is important mostly for the photon en- 
ergies hv < 1.0 keV. This indicates that a larger fraction of X-ray 
flux is affected/attenuated for the blackbody radiation curve with 
the smaller temperature (Tx — 2 x 10 K) since its peak is located 
below 1.0 keV. Consequently, the photoionization rate of He I is 
less affected for the higher Tx model (Model B6); hence, the line 
profile for He I A10830 is less affected for the higher Tx. 

From this test, it is clear that the assumption of no X-ray at- 
tenuation could significantly overestimate the line emission and ab- 
sorption strengths in He I A10830. In this test, we have considered 
the models with the X-ray flux attenuation which is calculated from 
the opacity of the models with no X-ray attenuation; hence, they are 
not completely self-consistent. This method may be acceptable for 
obtaining rough estimates of the effect of the X-ray attenuation on 
the line profiles. However, to implement the X-ray attenuation ef- 
fect properly in our radiative transfer model (KU1 1), the radiation 
field (including the X-ray radiation) and opacity must be solved si- 
multaneously. This requires an iteration between the opacity and 
the radiation field calculations until they are mutually consistent. 
Unfortunately, our radiative transfer code (KU11) is not designed 
for this type of iterations, and we have to assume that the gas is op- 
tically thin for the underlying continuum source. This is a limitation 
our radiative transfer model and is an important caveat. Therefore, 
we warn readers that the emission and absorption strengths of He I 
lines presented in this work could be overestimated significantly 
due to the assumption of no X-ray attenuation in our model. 

However, in reality, the magnetospheric accretion funnel, 
which is a significant contributor to the X-ray optical depth, may 
not be axisymmetric, and it could occur in two or more funnel 
streams, as shown in 3D MHD simulations (e.g. Romanova et al. 
2003 ). This flow geometry may contain large gaps between the fun- 
nel flows, and it could allow a large fraction of the X-ray produced 
near the stellar surface to reach the conical-shell regions without a 
significant attenuation, and ionize He I in the wind. 



5.3 Dependency on inclination angle and additional lines 

We now examine the dependency of the model profiles on the in- 
clination angles (i) and study the general characteristics by us- 
ing the set of model parameters used for Model B2 (see Ta- 
ble^, i.e. T max = 1.0 x 10 4 K, i x = 4x lO^ergs" 1 and 
Tx = 2 x 10 6 K. In addition to He I A10830 and Ha, samples 
of He I A5876, H/3, Pa/3 and Pa7 line profiles are also computed. 
Since the accretion disc region (8 <; 65° as shown in Fig. BJ is 
excluded in the opacity and emissivity calculations, the line pro- 
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Figure 5. Dependency of He I A 10830 and Ha profiles on the X-ray lumi- 
nosity Lx ■ The MHD simulation, Model B, is used for all the models shown 
here. The profiles are computed at the fixed inclination angle i = 30° 
and fixed X-ray emitting plasma temperature Tx = 2 X 10 6 K (assumed 
isothermal). The values of Lx used are 4 X 10 29 (solid), 4 X 10 30 (dash) 
and 2 X 10 31 erg s^ 1 (dash-dot), which corresponds to Models B4, B2 and 
B5 in Table pi respectively. The X-ray luminosity between 4 X 10 29 and 
4 X f 31 erg s _1 produces reasonable line strengths in He I A10830 when 
compared to observations (e.g. ED06). The model Ho profiles are not af- 
fected by the change in the values of Lx. le - three line profiles are almost 
identical to each other. 
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Figure 7. Demonstration of the effect of the X-ray opacity on the formation 
of He I A10830. The He I A10830 line profiles from Models B2 (T x = 
2 X 10 6 K) and B6 (T x = 5 X 10 6 K) (solid) are compared with the 
same models computed with the input X-ray continuum fluxes corrected 
for the attenuation by the continuum optical depth in X-ray (dash). Note 
that the original Models B2 and B6 have the same X-ray luminosity, Lx = 
4 X 10 3 " erg s — x (Table|2l. For a relatively low X-ray temperature (Tx = 
2 X i0 6 K) (left panel), the reduction of the line emission and absorption 
strengths due to the attenuation of the X-ray flux is significant. The effect 
of the X-ray attenuation becomes smaller for the model with a higher X-ray 
temperature (Lx = 5 X 10 6 K) (right panel). 
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Figure 6. Dependency of He I A 10830 and Ha profiles on the X-ray temper- 
ature Tx. The MHD simulation, Model B, is used for all the models shown 
here. The profiles are computed at the fixed inclination angle i = 30° and 
with the fixed X-ray luminosity Lx = 4 X 10 3() erg s — 1 . The values of Tx 
used are 2 X 10 6 (solid), 5 X 10 6 (dash) and 2 X 10 7 K (dash-dot), which 
corresponds to Models B2, B6 and B7 in Table[2] respectively. The X-ray 
temperature between 2 X 10 6 and 2 X 10 7 ergs _1 produces reasonable 
line strengths in He I A10830 when compared to observations (e.g. ED06). 
The model Ha profiles are not affected by the change in the values of Tx, 
i.e. three line profiles are almost identical to each other. 



files presented hereafter are restricted to i < 65°. This is done 
also to avoid the complication that arises from the occultation of 
the line emission by the extended disc which we do not treat self- 
consistently. We choose four different inclination angles, equally 
spaced in cosi between i — 0° and 65°, to present the profiles 
of equal probability of occurrence between the inclination angle 
range, i.e. they are approximately i = 22°, 38°, 50° and 60°. The 
results are summarized in Fig. [8] In general, the line strengths of 
the model profiles are comparable to those found in observations 
(e.g. |Reipurth et al.|[T996; Edwards et al.||1994||Muzerolle et al.| 



|1998b||Alencar & Basri|2000[|Folha & Emerson|2001| ED06), ex 
cept for He I A5876 which are much weaker than a typical value 
from observations (see e.g. Beristain, Edwards, & Kwan 2001 ). In 
the following, the characteristics of each line are examined more 
closely. 

He I X10830. The conical-shell wind geometry with the half- 
opening angle ~ 35° (Section |4.2.2[ > does not favour the wind 
absorption at a low inclination (e.g. i = 22°). The blueshifted 



absorption, caused by the conical-shell wind, starts to appear at 
i « 30°. The position of the blueshifted absorption shifts from 
V ~ 120 km s _1 (at i « 38°) towards the line centre as the in- 
clination becomes higher (up to i « 60°). This shift is caused by 
the change in the line-of-sight velocity of the conical-shell wind 
toward the observer as the inclination changes, The redshifted ab- 
sorption, caused by the accretion funnel (c.f. Fig. RJ, is present 
for all i. The extent of the redshifted absorption changes between 
~ 200 to ~ 250 km s _1 depending on i, which is similar to the 
flow speed of the gas in the accretion funnel near the stellar surface 
(c.f. Section |4.2.2| >. At higher inclination angles (i <; 50°), the line 
emission almost disappears, and the profile is shaped by the wind 
and funnel absorption features. 

He I \5876. The line peak flux slightly decreases and the red- 
shifted absorption component becomes stronger as the inclination 
increases. No blueshifted wind absorption component is seen in this 
line for all i. This is in agreement with the observation of Beristain 
|et al.| ( (2001[ >. The redshifted absorption components, on the other 
hand, are present for all i. However, the observations of this line, 
from 31 CTTSs ( |Beristain et~aLl[200T] |, show that the redshifted 
absorption is rather rare (3 out 31). The line peak fluxes in our 
models are between ~ 1.1 and ~ 1.2 (normalised to continuum) 
which are much smaller than those found in the observations (typi- 
cally > 1.2, e.g. |Alencar & Basri|2000||Beristain et al. 2001 1. This 
line is thought to arise from the post shock region of the accretion 
funnel near the stellar surface (e.g. Beristain et al. 2001 , Edwards 



|et al.|2003|rOrlando et al.|2010||Kwan & Fischer|201 \) . This indi- 
cates that a higher temperature and density gas may be needed in 
order to match the observations which would require more formal 
treatment of the shock regions in the accretion funnels (e.g. Lamzin| 



[19981 [Koidbba et al.|2008"l|Sacco et al.l2008"l|Orlando et al.|2010> " 
Ha. Unlike the helium lines, no clear redshifted absorption 
component is seen in this line. However, the profile at the low in- 
clination (i — 22°) shows a slight blue asymmetry, i.e. the flux 
distribution is slightly shifted toward blue, which may be caused 
by a small amount of absorption in the red wing although it is not 
seen as a discrete absorption component. At a low inclination an- 
gle (i ss 22°), no clear blueshifted wind absorption component 
is seen in the model profiles. A wind absorption becomes visible 
for i <; 38°, and its strength increases as i increases. In general, 
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Figure 8. Summary of sample model hydrogen and helium line profiles computed at four different inclination angles: i = 22°, 38°, 50° and 60° (equally 
spaced in cos i between i = 0° and 65°) from the top to bottom rows, respectively. From the left to right columns, sets of He I A10830, He I A5876, Ha, H/3, 
Pa/3 and Pa7 profiles are shown respectively. The X-ray luminosity Lx = 4 X 10 30 erg s _1 and temperature Tx = 2 X 10 6 K (Model B2 in Tablepl are 
used for all the models shown here. The fluxes (F n ) are normalised to the local continuum. 



the models are in good agreement with observations (e.g.|Reipurth 



et al. 1996; Edward s~eTaT1[T994l |Muzerolle et al"l|1998b| |Alen 



car & Basri 2000). The model profiles here belong to Type I and 



Type III B (which constitute about 60 per cent of the observed Ha 
profiles for CTTSs in Reipurth et al. 1996) according to the mor- 
phological classification of Reipurth et al. ( 1996). 

Hj3. The redshifted absorption component is present in the 
profiles at all i. The extent of the redshift absorption changes 
between ~ 200 to ~ 250 km s _1 depending on i. While no 
blueshifted absorption component is seen at the lower inclination 
angles (e.g. i = 22°), a strong blueshifted wind absorption com- 
ponent is present at higher inclinations (i > 38°). The position 
of the wind absorption component shifts toward the line centre 
as i further increases. At i — 60°, the wind absorption compo- 
nent is located near v = Okms -1 , and the resulting line has a 
double peaked profile shape. At a very high inclination angle, the 
conical-shell wind is moving almost perpendicular to the line of 
sight to an observer; hence, the wind absorption is located near 
the line centre. The double-peak profile resembles that from a ro- 
tating disc; however, the H/3 emission here mainly arises from the 
magnetospheric accretion funnel whose rotational velocity is rather 
slow (~ 20 km s^ 1 ). The extent of the emission in the blue wing 
becomes larger as i increases. This is consistent with the trends 
found in the radiative transfer models of Hart mann et al.| (1994) 
and Muzerolle et al. (2001 ) who modelled the line emission from 
the axisymmetric magnetospheric accretion funnel flows. The line 
profile shapes and the strengths found here are comparable to those 
found in the observations (e.g. Edwards et al. 1994; |Alencar &] 
|Basrip 000). However, the extent of the blue wing in the model 
profiles (~ —130 km s _1 to ~ —200 km s _1 ) tends to be notably 
smaller than a typical value found in observations (~ 300 km s _1 ). 

Pa/3. The inclination angle dependency of the line profiles are 
somewhat similar to that of H/3 and Pa7. The redshifted absorption 
component is visible at all i, and its depth slightly decreases as i in- 
creases. The extent of the line emission in the blue wing increases 
from ~ —130 km s _1 to ~ —200 km s _1 as i increases from 22° 



to 60°. At i = 22° and 38°, the models show a classic inverse P- 
Cygin profile, which is classified as Type IV R profiles according to 
Reipurth et al. ( 1996) and Folha & Emerson (2001 1. This consists 
of about 34 per cent (13 out of 37) of the objects observed by |Folha| 
& Emerson (2001). The blueshifted wind absorptions are seen at 
the higher inclinations (i — 50° and 60°), and their positions shifts 
towards the line centre as i increases. Similar to the H/3 case, the 
wind absorption component is located near v ~ km s _1 for the 
model with i — 60°, and the line exhibits a double peaked shape. 
On contrary, the blueshifted wind component is extremely rare (1 
out of 37 objects) in the observed Pa/3 profiles of Folha & Emerson 
(2001). This indicates that the Pa/3 line opacity in the conical-shell 
wind is too high possibly because the wind density and/or temper- 
ature are too high in our model. 



Paj. The dependency on i is very similar to that of Pa/3 and 
H/3. The line shows the redshifted absorption component at all i. 
The frequency of the redshifted absorption component is much 
lower (24 per cent) in the observation of ED06. A relatively strong 
blueshifted wind absorption appears in the higher inclination an- 
gles, i.e. i > 50°. The peak fluxes and the profiles shapes of the 
models are similar to those found in the observations of 38 CTTSs 
obtained by ED06, except for the blueshifted absorption feature in 
the models with i > 50°. No wind absorption component is seen 
in the observed Pa7 profiles of ED06. As in the Pa/3 models, this 
indicates that the line opacity in the conical-shell wind is too high 
possibly because the wind density and/or temperature are too high 
in our model. The extent of the blue wing emission increases from 
~ — 130kms~ to ~ — 200 km s -1 as i increases from 22° to 
60°. This behaviour, again, is consistent with the line profile mod- 
els of |Hartmann et al.| l |1994) and |Muzerolle et al.| ( |200l) . 
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6 DISCUSSION 

6.1 Contribution to line emissions from the conical-shell 
wind 

The model line profiles presented so far are the emergent profiles 
which an observer at a given viewing angle would observe. This 
emission consists of the emission from the inflowing gas (the ac- 
cretion funnel), the outflowing gas (the wind), and the continuum 
emission. Although the red and blueshifted absorption components 
in the model line profiles could provide us some information on the 
physical conditions and locations of the line 'absorption', the loca- 
tions for the line 'emission' could not be easily interpreted from the 
total emergent line profiles. For this reason, we now examine the 
relative contributions of the accretion funnel and the conical-shell 
wind to the emergent profiles by comparing the line profiles com- 
puted with (1) the accretion funnel part only, (2) the conical-shell 
wind only, and (3) the combination of (1) and (2). The physical pa- 
rameters used for this demonstration are the same as in Model B2 
(see Tablepl. All the line profiles here are computed at the inclina- 
tion angle i = 30°. The results are placed in Fig. [9] 

The difference in the continuum levels for each line profile 
is very small, i.e. the continuum emission from the wind and the 
accretion funnel is negligible compared to the stellar continuum; 
therefore, the line profiles presented in Fig. [5] use the normalized 
flux for the vertical scales. The figure shows that, with this partic- 
ular set of model parameters, the emission in He I A10830 mainly 
arises from the accretion funnel (where the photoionization radia- 
tion is very strong); however, the emission from the wind is non- 
negligible, and contributing to the line centre flux. On the other 
hand, the wind emission in He I A5876 is negligibly small, and the 
line emission is essentially from the accretion funnel. Unlike the 
helium lines, the wind emission in the hydrogen lines is comparable 
to that of the accretion funnel for H/3 and Pa7, and is significantly 
dominating that of the accretion funnel for Ha and Pa/?. However, 
the readers must be aware that the wind emission in hydrogen lines 
are very sensitive to the assumed wind temperature which is essen- 
tially controlled by our model parameter T max , as demonstrated in 
Section |5T| for Ha. The wind temperature in these models is ap- 
proximately isothermal with T ~ T max = 10 4 K; however, if a 
slightly lower wind temperature, e.g. T max = 8 x 10 3 K is used, 
the contribution from the wind emission will decrease significantly. 

Despite the uncertainty in the wind temperature, it is inter- 
esting to note that the extent of the blue wing in some hydrogen 
lines, especially Pa/3 and Pa7, becomes notably (~ 50kms _1 ) 
larger due to the wind emission. However, the amount of the blue 
wing extention would depend on the relative flow speed of the gas 
in the accretion funnel and the conical-shell wind. Since the wind 
emission in the hydrogen lines presented here is significant, the fi- 
nal emergent line profiles (computed with both wind and accretion 
funnel) appear broader than those from the accretion funnel or mag- 
netosphere alone. However, this also depends on the assumed wind 
temperature. 



6.2 Comparisons of He I A10830 profiles with observations 

The sensitivity of He I A10830 to the innermost winds of CTTSs 
and its usefulness for probing the physical conditions of the winds 
have been demonstrated by e.g. ED06, |Kwan et alT|j2007) , |Kwan| 
|& Fischer] ( |2011| > and KU11. Here, we examine whether the 
conical-shell wind model (Model B) can account for the types of 
He I A10830 profiles seen in the observations. For this purpose, we 
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Figure 9. Comparison of the model lines profiles (He I A10830, 
He I A5876, Ha, H/3, Pa/3 and Pa7 ) computed with different flow elements 
in the MHD simulation (Model B). For each line, the profiles are computed 
with (1) the accretion funnel part only (dash), (2) the conical-shell wind 
part only (dash-dot), and (3) all the flow components - the accretion funnel 
and the wind included (solid). The physical parameters used here are the 
same as in Model B2 (Tablep), and the inclination angle is i = 30°. Note 
that the flux contribution from the conical-shell wind in the hydrogen lines 
is sensitive to the assumed wind temperature which is essentially controlled 
by the model parameter T max (see Section [5~[) . The contribution from the 
wind would decrease significantly when T max is decreased. 



have run a small set of line profile models for a various combina- 
tions of Lx, Tx and i while keeping other parameters fixed as in 
Section [531 We then examined if there is a resemblance between 
any of our model profiles to a set of He I A10830 line observations 
in ED06. In the following, we present a several example cases in 
which the model profile morphology is similar to that of observa- 
tions. 

Fig. [To] shows qualitative comparisons of the model 
He I A 10830 profiles with the observation of ED06. The corre- 
sponding model parameters used are summarised in Table [3] The 
models shown here are not the strict fits to the observations, but 
rather simple comparisons of the profile morphology since we did 
not adjust the underlying stellar parameters (stellar masses, radii, 
luminosity and so on) which could be tailored to individual objects. 
The figure shows that the line profiles predicted from the conical- 
shell wind model (Model B) broadly agree with the observations. 
The extent of the redshifted absorption, formed in the accretion fun- 
nel flow (cf. Fig. Mb, in the model line profiles is ~ 250km s _1 , 
which is comparable to those in the observed profiles (~ 200 to 
~ 300 km s _1 ). The conical-shell wind produces a relatively nar- 
row blueshifted absorption component in the profile in various de- 
grees depending on the model parameters, Lx, Tx and i. As we 
demonstrated in the previous section (Figs.[5][6]and[8]in SectionBl, 
the wind absorption depths and the emission strengths are sensitive 
to these parameters. 

A relatively strong line centre emission (as seen in the case for 
CY Tau) can be obtained with a relatively high X-ray luminosity 
(Lx = 8 x 10 30 ergs~ 1 ) and a relatively low inclination angle 



© 0000 RAS, MNRAS 000, 000-000 



14 Kurosawa & Romanova 



Table 3. Summary of models used for comparison of He I A 10830 with 
observations in Fig.fTol 



Object Lx Tx inclination 

(lO^ergs" 1 ) (10 6 K) ••• 



Model 



Observation 



CYTau 

CITau 

HKTau 

UZ Tau W 



8 

2 

2 

0,1 



2 30° 

2 36° 

5 32° 

2 30° 



(i — 30°). On the hand, a relatively weak line centre emission 
and wind absorption (as seen in CI Tau, UZ Tau W and HK Tau) 
can be obtained with a relatively low X-ray luminosity(Lx = 0.4— 
2 x 10 30 ergs -1 ) and a relatively low inclination angle (i — 30°- 
36°). A notable difference between the models and observation is 
in the extent of the blueshifted absorption component in UZ Tau W. 
In both cases, the observations suggest a presence of high velocity 
outflows (up to ~ 300 km s _1 ) which cannot be explained by the 
conical-shell wind in our MHD model because its maximum speed 
is only ~ 125 km s _1 (Fig.pk. The high velocity component seen 
in the observations are possibly attributed to a bipolar 'stellar wind' 
which are not implemented in our model. 

As demonstrated above, the conical-shell wind model 
(Model B) are capable of reproducing the relatively narrow and 
low-velocity (v ~ —130 to —200 km s _1 ) blueshifted absorption 
component in He I A10830 profiles. However, the model does not 
show the profiles with a deep and wide blueshifted wind absorp- 
tion component that extends to v ~ —300 to —400 km s _1 , which 
are observed in about 40 per cent of the samples in ED09. These 
P-Cygni like line profiles are likely caused by stellar winds that 
emerge in the polar directions of CTTSs, as suggested by ED06, 
|Kwan et al.| |2tJ07) , |Kwan & Fischer) pOTT) and KU1 1 . The in- 
ability of our model for reproducing the P-Cygni like profile was 
expected since our MHD model does not include a stellar wind. 
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Figure 10. Qualitative comparison of the model He I A 10830 profiles with 
the observation of ED06. Example model profiles (left column) are selected 
based on their morphological similarities with the observations (right col- 
umn). For the observed profiles, the names of the objects are indicated on 
the upper-right corner of each panel. The corresponding model parameters 
used are summarised in Table [3] The He I A 10830 profiles predicted from 
the conical-shell wind model broadly agree with the subset of the observa- 
tions, but the model does not show a deep and wide blueshifted (P-Cygni 
like) profiles (not shown here) found in some CTTSs (see ED06). 



7 CONCLUSIONS 

We have presented the line profile models of hydrogens and he- 
lium based on the axisymmetric MHD simulation of the outflow 
produced at the interface of the magnetosphere of CTTSs and their 
accretion discs. The so-called conical-shell wind (RO09) is formed 
when the stellar dipole magnetic field is compressed by the accre- 
tion disc into the X-wind like configuration. We have extended 
the previous conical-shell wind simulations of RO09 to include 
a more well defined magnetospheric accretion funnel flow which 
is essential for modelling the optical and near-infrared hydrogen 
and helium lines of CTTSs, as many of the lines show clear signs 
of inflow (seen as the inverse P-Cygni profiles). The formation 
of the well defined funnel flow with its outer radius correspond- 
ing to ~ 2 Rq (Model B in Section |4.2} is obtained by lowering 
the disc density or correspondingly the mass-accretion of the disc 
from the previously found conical-shell wind solution of RO09 in 
which the magnetosphere was compressed all the way to the stel- 
lar surface by the accretion disc. The mass-accretion and mass-loss 
rates (in the conical-shell wind) found in the model (Model B) are 
4.1 x 10" 8 Mq yr -1 and 6.9 x 10" 9 M Q yr -1 , respectively (Ta- 
ble [TJ. The maximum velocities in the conical-shell wind and the 
funnel flows are ~ 125 km s _1 (Fig. |2p and ~ 250 km s _1 , re- 
spectively. These values are all comparable to the values found in 



observations (e.g. Hartigan et al. 1995; Gullbring et al. 1998 Cal- 
|vet et al.||2050l [Alencar & Basri||20001 |Folha & EmersonpOOl] 
ED06). In the following, we summarise our main findings from the 
line profile models based on this MHD simulation. 

We find that He I A10830 is relatively sensitive to changes 
in the X-ray luminosity (Lx) and X-ray temperature (Tx) (Sec- 
tion ^. 2} . Both emission and absorption components of the line be- 
come stronger as Lx increases (Fig.Bl. For a fixed Lx value, the 
emission and absorption components of the line become weaker as 
Tx increases (Fig.[6l because the peak of the thermal radiation (as- 
sumed as a single temperature gas) shifts to a higher energy and 
the amount of the soft X-ray emission decreases. The softer X-ray 
photons photoionize more He I since its photoionization cross sec- 
tion is larger towards the EUV frequency. The range of Lx between 
4 x 10 29 and 2 x 10 31 erg s _1 produces reasonable line strengths in 
He I A10830. Similarly, T x between 2 x 10 6 and 2 x 10 7 K (with 
L x = 4 x 10 30 ergs _1 fixed), also produced He I A10830 with 
its line strength comparable to observations. On the other hand, yet 
another wind sensitive line, Ha, shows no significant change in its 
line strength when Lx and Tx are varied in the same ranges used 
for He I A10830 because the photoionization rates (in the X-ray 
range for H I) are small compared to the collisional rates. 

A rich diversity of line profile morphology is found in our 
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sample hydrogen and helium model profiles (Fig.[8]in Section [53] l. 
Many of the model profiles are very similar to those found in the 
observations (e.g. Reipurth et al. 1996; Edwards et al. 1994; Muze- 
|rolle et al.|J998bl|Alencar & Basri|2000||Folha & Ernerson|2001| 
ED06). No wind absorption is seen at a very low inclination an- 
gle (e.g. i < 22°) in any of the lines considered here. In addi- 
tion, the wind absorption component seen in the higher inclina- 
tion angles are rather narrow with the corresponding wind speed 
v Ji 120 km s" 1 . These suggest that the blueshifted absorption 
components seen the model profiles are caused by the conical-shell 
wind and no significant absorption occurs (at least with the model 
parameters used here) in the low-density polar wind. The main 
shortcomings of our models are in the weakness of He I A5876 
emission, and in the presence of the strong wind absorption com- 
ponent in Pa/3 and Pa7 which are not seen in the observations. 
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We have examined the relative contributions of the conical- 
shell wind and the accretion funnel to the hydrogen and helium 
line emissions (Section |6. 1) . We find that the conical-shell wind 
can significantly contribute to the line emissions in Ha, H/3, Pa/? 
and Pa7, provided that the temperature in the wind is high enough 
(e.g. ~ 10 4 K). On the other hand, the helium lines (He I A10830 
and He I A5876) have less contribution from the wind emission. 
The wind contribution to the hydrogen emission decreases signif- 
icantly when a lower wind temperature is adopted in the profile 
calculations. A better constraint of wind temperatures could be ob- 
tained from the comparisons of the line ratios of e.g. the Balmer 
lines from observations and that from line profile models. 

Using the various combinations of our main model parame- 
ters (Lx, Tx and i), we have demonstrated that the conical-shell 
wind model is capable of producing a subset of the line profile 
morphology for He I A10830 found in observations (Fig,|10[>, The 
model is able to reproduce a relatively narrow and low-velocity 
blueshifted absorption component in He I A10830 profiles, found in 
about 30 per cent of the sample in ED09, provided that the inclina- 
tion angle is favourable (~ 30° to ~ 50° for Model B, but the range 
will depend on the opening angle of the conical-shell wind). How- 
ever, the model does not reproduce the deep and wide blueshifted 
wind absorption components (the P-Cygni like line profile) found 
in about 40 per cent of the sample in ED09. The earlier studies by 
|Kwan et ai] ( |2007^ and KU1 1 demonstrated that the P-Cygni like 
profiles seen in He I A10830 of CTTSs are most likely produced by 
the absorption by the stellar wind that blows in the polar directions. 
The inability of our model for reproducing the P-Cygni like pro- 
file was expected since our MHD model did not include the stellar 
wind component. 

Finally, in a future study, we plan to implement the stellar 
wind in our MHD simulations, and study its effect on the line for- 
mations, especially on the P-Cygni like profile of He I A 10830. We 
also plan to investigate the line formation problems in the higher 
mass-accretion regimes which are somewhat similar to Model A 
(Fig.fTlin SectionBl, and those in Konigl et al. (201 1 1 and Lii et al. 
(2012). These models, in which the compression of the magneto- 
sphere toward the stellar surface is strong and the outflow due to 
the magnetic pressure gradient force is very strong, may be appli- 
cable to the objects during the outbursts periods such as EXors and 
FUORs. The time-series line profiles of the outburst can be mod- 
elled by using multiple time slices (rather than a single time slice 
used in this work) of MHD simulations, to study the time evolution 
of observed spectra. 
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